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)^', We study a general field theory of a scalar field coupled to gravitation through a quadratic 

Gauss-Bonnet term C,{4>)Rgb- We show that, under mild assumptions about the function 
C,{4>), the classical solutions in a spatially fiat FRW background include singularity - free 
solutions. 
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Recently it was found that the loop corrected superstring effective action in a FRW 
background in the presence of the dilaton and modulus fields admits a particularly interesting 
class of singularity-free solutions with flat initial asymptoticsf^'. Singularity-free solutions have 
been considered by other authors as welP'''^^ and their existence, although linked to the 
presence of quadratic curvature terms, depends on the form of the gravitational coupling to 
scalar fields. The purpose of this letter is to study the existence of singularity- free solutions 
and formulate the necessary conditions that should be imposed on the quadratic gravitational 
coupling. Although motivated by the concrete case of the superstring the analysis at hand is 
quite general and is carried out in an effective field theory framework that allows a possible 
alternative, at present non-existing, quantum gravitational underlying theory. 

Consider a scalar field coupled quadratically to gravity through the Langrange density 
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L = -R+-iD<pf--a<P)Rl^ 



where Rqb ~ {^hukx)'^ — '^{RjivY + R^ is the standard Gauss-Bonnet integrand. The function 
^(0) is a general smooth function of the scalar field which will be further constrained shortly. 
5 is just a coupling parameter. Note that we have chosen a Gauss-Bonnet quadratic coupling 
in accordance with general unitarity arguments as well as the concrete superstring case^. 
In a standard RWF background metric g^y = (1, — e^'^^jj) we obtain the classical equations 
of motion 

0" + 3u<P + pei<P)^\^" + dj^) = (2) 

3^2 _ 1^2 _ 3^^.(^)^^3 ^ (3) 

In terms of x = and z = Co considered as functions of 0, equation (3) can be solved as an 
algebraic equation and give 



while (2) becomes 



x = z[-^6^'z'±J6 + {pez^y] (4) 



z'^^ = --iA/B) (5) 

ad) X 



where 

A = 16/ + IQz^x^ + -x^ - hCz'^x" (6) 
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B ^ 16/ - l^/a;2 + ^x^ (7) 
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Let us assume now that ^(0) is a smooth function that possesses one minimum at some 
point 00 • For example ^ could be 0^,0*^,... . 

We shall show that under the assumption that .^(0) grows faster than 0^ at ^ ±oo 
the system of equations of motion (2), (3), admits non-singular solutions for the case 5 > 0. 
Let us outline here the basic steps of our proof. We first show that all singularities occur at 
\z\ -^ oo. Expanding asymptotically our equations at this region and analyzing the various 
possibilities we show that only one singular solution exists oj ~ logt, ~ 0^ + t"^ / {6 (,'{(/) s))- 
Then we show that the points z = or x = cannot be continuously approached for finite 
and thus every solution is characterized by a definite sign of z and x. Choosing z > 
(expanding universe) and a: > we show that the singularity point 0s is related to the ^ 
minimum 0o, and we always have 0^ > 0o, for 6 > 0. Furthermore, we show that a singular 
solution involves only points > 0^ and thus any solution that includes a point < 0o 
is non-singular. Since such a point can be always found this is an existence proof of the 
singularity free solutions. Finally, in order to make the results of our analysis concrete, we 
present the numerical solutions of the case C, = 0^/2, 6 = 1. 

A "singular solution" is one for which \z\ ^ oo or |0| ^ oo at some finite time. Note 
that we could have z ^ oo for — ^ oo and still have a singularity at a finite time. We shall 
restrict the coupling function ,^(0) to increase at least as fast as 0^ at infinity. Then, we can 
show that x{(f)) cannot have a singularity while z is maintained at a finite value. In order to 
prove it we first solve equation (4) for large x and obtain x ~ —3(5.^ z^ . Obviously x has a 
singularity while z is finite, only if 6^ (0) has one. Since ^(0) is supposed to be a smooth 
function this can occur only at ^ ±oo . Next, we go to equation (5) which takes the form 

/ . 1 ... 6 
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Sot, z 5 

We can proceed considering three separate cases. Namely SC, z^ ^O(l) , 5S, / ~0(1) 
and 5^ z^ ^O(l) while always b^ z^ ^O(l) . Solving (8) in the first case leads to the 



contradiction x ~ {6C, )"^ -^ 0. The second case can be treated by setting 5^ z^ = h 
and getting from (8) the solution z ~ 02(i~''), % ~ 02~f , S^ ~ 0^"^'' for -^ ±oo with 
K = (| + ;^)^"'^- In order to have the assumed behavior of z^x and 5^ we have to constrain 
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1 < K < 3. Solving in terms of time we obtain ~ (t + ...) ("-i) . Thus, the point — * ±oo 
is not approached at a finite time and we do not have a true singularity in this case. Finally, 
in the last case 5^ z^ <^0(1) the equation (8) becomes z ~ —^ which is satisfied by 
z^b^ ~ |0 for -^ ±oo . Expressing (8) solely in terms of z gives z'^ ~ 0, b^ ~const. 
which contradicts our assumptions. 

Let us now proceed by deriving the asymptotic form of the singular solutions. We can 
start from equation (4) and solve it around the singular point 0s for which -2(0^) -^ oo . 
Three different possibilities arise in principle for z near the singularity, namely b^ z^ ~0(1), 
b^ z^ ^O(l) and b^ z^ ^O(l) . They can be analyzed as follows: 

1) b^ z^ ~0(1). In this case we set |(5^ z^ = a near 0^ and obtain from (4) x = \z with 
a = I — |. Equation (5) gives 

/ -(l + A' + t-^) 

A(l-f + t) '"' 

or ^ = —^ ■ For b^" z"^ = b ~0(1) we get z ~ exp(— i^) which is singular at ^ ±oo 
and leads to b^ ~ 6exp(— ^^) in contradiction to our assumptions for the behavior of ,^(0) 
at infinity. For b^"z'^ ^O(l) we get the solution z^"^ ~ -^ + c which requires c=0 and 
| — 1 = 9(1 — ^ + ^)/2A^. The last relation is impossible since there are no real solutions 
for A. Finally, the case bC,"z^ ^O(l) is just like the first case. 

2) b^'z"^ <^0(1). In this case from equation (4) we obtain x ~ ±z^/Q . Equation (5) 
becomes 

(~T3(l-^)/v^ (10) 

If b^ z"^ = b ~0(1) we obtain z ~ exp(=l=^/|(l — |)0) which corresponds to a singularity 



only for -^ ±oo . This implies b^ ~ exp(±2w|(l — |)0) which is not in agreement with 
our assumptions for ,^(±oo). Similarly for the case b^ z^ ^O(l). Finally, if b^ z^ ^O(l) , 
(10) gives z^'^ ~ ±w|5.^ which is incompatible with b^ z^ ^O(l). 



3)6C, z^ ^O(l). In this case we get two different solutions for x, namely X- ~ —3(5^ z^ 
and x^ ~ Sr- ■ 

For X = X-^ neccessarilly |a;| ^ \z\ and equation (5) becomes 

' ^ ^^^ (11) 



which is equivalent to (j^b^ z^^ ~ 1 + ^b^' z^ . If b^' z^ = b ~0(1) we obtain ^b^' z'^ ~ 
(1 + 1^)0 + ... which implies that the singularity appears for infinite . Then, we get 



b^ ~ 0^^/(2+V) and z'^ ~ 0(i^~)/(i+to)j with 6 < | in order to have |z| ^ oo . The constraint 
6 < I leads to a function bC, that does not increase at infinity as fast as 0^, in contradiction 
to our assumptions. If b^ z"^ ^O(l), we obtain ^bC, z^ ^ (p which also requires -^ ±oo . 
Then, we are led to b^ 4>/bC,' ^O(l) which is not acceptable. Finally, for b^ z^ ^O(l) we 
get from (11) z ~ [b^ )^5 which is not consistent with z^b^ ^ 1 at z ^ oo. Thus, the case 
of X- is impossible. 

For a; = a;+ ~ -^ we must have \xl z\ <^ 1 . Then, equation (5) becomes 



SXC" 



If b^' jz'^ib^f <0(1), we obtain 

"^ ^ bm - m<ps) ^^^^ 

which is a singular solution at some finite point 0s • If b^ /{b^ )^ ^ z^, we can get from 
equation (12) z"^ ~ ~^/{{^^ )^ + •••)) which is impossible as 2; ^ oo . The only remaining 
case ,^|,,4 ~ 0{z^) can be studied by putting b^" = bz^{b^'y. Then, from (12) we get 
2;-2 ^ ^i _ ^-^^(. _^ ^^^^-^y Solving for b^' we also get {b^')^^ ~ const. + gj, 3 , which 

is equivalent to ^ r' iJ, ~ 1 • This is however, for any non-zero b, contradictory to our 

(1^-3") 

assumptions. 

Thus, assuming that .^(0) is a smooth function that grows as fast as 0^ at both infinities 
-^ ±00, we have found only one singular solution. Near the singular point 0s {C, (0s) 7^ 0) 
it behaves as 
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z r^ {Sm - 5a<Ps))-- (14a) 

X - j^ (146) 

In terms of the time, (14) becomes 

<P^ <Ps + -FTTT-T (15a) 

0^ (0s) 

cj~lnt (156) 

where we have shifted the time variable so that the singular point (ps corresponds to t=0. 
Note that while z ^' oo, x goes to zero as long as 6C, (0^) 7^ . 

The limit x ^ of the system (2), (3) deserves some more attention. Indeed, as can 
be seen in equation (4), when a; ^ for finite the only possibility for z is to go to zero 
too while X ~ ±^/6z . Substituting the last relation in equation (5), we obtain z j z r^ 
"Fyl or z ~ exp=l=w|0 which does not go to zero for finite in contradiction with our 
hypothesis. Inversely, when 2; — * for finite the only possibility read off from equation 
(4) is a; ~ ±v^2;. Again, substituting that in equation (5) we reach the same contradiction. 
The only conclusion we can draw is that the points x = or 2; = cannot be reached 
analytically. Therefore, the signs of x and z must be independently conserved. Then, the 
following statement is true for the singular solution we found: Every connected piece of 
the solution will be characterized by a fixed sign for each of x and z . We have assumed 
single- valuedness of -2(0). If in a connected piece of the solution z(0) were not single- valued 
then at a non-singular point 0i we would have ^ = at 2:1 = z{(\)\) . This is equivalent 
to xB j Az\ = and it can happen only when either a; -^ or A ^ 00 in contradiction to 
Zi < 00. 

Consider now a singular solution with the singularity occurring at a point 0s . Each 
connected piece of the solution will be characterized by a given sign of x and z . Take the 
piece with x > 0, z > . Then, 

x = ^--^ (16) 



implies that SC, > and therefore, ii S > 



>0O 



where 0o is the minimum of ^(0) . Because of z ~ —^z^6C, < and the single-valuedness of 
z we come to the conclusion that these solutions cover at best the full quarter-plane z > 0, 
> 00 but they do not contain any points of the z > 0, (f) < (po region. We are now in a 
position to state the following "theorem" : 

For 6 > 0, any solution containing a point 0* < 0o such that x(0*) > 0, z(0*) > must 
be a non-singular solution. 

Since the solution is controlled by one first order differential equation, namely (5), while 
(4) is just a soluble algebraic equation, it depends only on z, x at some initial point and we 
can always satisfy the requirements of the theorem starting with values z* > 0, x* > at 
some initial point 0* < 0o . 




Figure 1: The phase space diagram z(0) for the case ^(0) = 20^, 6 = 1 and 2; > 0, x > 0. 
The dashed lines correspond to singular solutions and the continuous lines to non-singular 
ones. Arrows indicate the direction of the time. 



The existence of the non-singular solution depends crucially on the fact that the singular 
solutions do not cross over to the (p < (po region. Notice that if we keep x > 0, z > and 
chose 6 < equation (16) demands (ps < 4>o ■ Then, since z < still, the whole upper plane 
is filled with singular solutions and there are no singularity-free solution in this case. This 
behavior is in agreement with the violation of the energy conditions required by singularity 
theorems ^^' which cannot occur in the 5 < case since in this case p + p and p + 3p are 
positive definite. 

In order to make our analysis concrete we have also studied numerically the case ^(0) = 
0^/2, 6 = 1. The phase diagram for z = z{(f)) > and x > is presented in Fig. 1. As 
expected the singular solutions can occur only for 0^ > 0o = and points with < 0o = 
belong to non-singular solutions. It can also be shown that in this specific case the non- 
singular solutions interpolate between a De-Sitter space universe at the remote past {t -^ 
— oo) and a slowly expanding one at the future {t -^ — oo). 

The general behavior of the non-singular solutions at late times {t —>■ oo) , assuming 
an expanding universe (2; > 0) , is that of slow expansion, while the behavior at early 
times(t — > —00) depends crucially on the particular choice of ^(0). In the case of the loop 
corrected superstring effective action '^^ which has served as a motivation for the present 
investigation we obtain at early times a flat-space. In contrast, in the case ^(0) = ^0^, 
analyzed above , we have obtained De-Sitter space. 
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^In the superstring effective action case [1], tlie modulus dependent coefficient ^{<j)) of tlie Rqs term, 
satisfies tfie requirements of our proof, since it has a unique minimum at </) = and it grows as ^ ~ el*^' at 
d> — > ±(X). 
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